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Abstract 

Continental flood basalts (CFB) are considered as potential C0 2 storage sites because of their high reactivity and 
abundant divalent metal ions that can potentially trap carbon for geological timescales. Moreover, laterally 
extensive CFB are found in many place in the world within reasonable distances from major C0 2 point emission 
sources. 

Based on the mineral and glass composition of the Columbia River Basalt (CRB) we estimated the potential of CFB 
to store C0 2 in secondary carbonates. We simulated the system using kinetic dependent dissolution of primary 
basalt-minerals (pyroxene, feldspar and glass) and the local equilibrium assumption for secondary phases 
(weathering products). The simulations were divided into closed-system batch simulations at a constant C0 2 
pressure of 100 bar with sensitivity studies of temperature and reactive surface area, an evaluation of the reactivity 
of H 2 0 in scC0 2 , and finally 1D reactive diffusion simulations giving reactivity at C0 2 pressures varying from 0 to 
100 bar. 

Although the uncertainty in reactive surface area and corresponding reaction rates are large, we have estimated the 
potential for C0 2 mineral storage and identified factors that control the maximum extent of carbonation. The 
simulations showed that formation of carbonates from basalt at 40 C may be limited to the formation of siderite 
and possibly FeMg carbonates. Calcium was largely consumed by zeolite and oxide instead of forming carbonates. 
At higher temperatures (60 - 100 C), magnesite is suggested to form together with siderite and ankerite. The 
maximum potential of C0 2 stored as solid carbonates, if C0 2 is supplied to the reactions unlimited, is shown to 
depend on the availability of pore space as the hydration and carbonation reactions increase the solid volume and 
clog the pore space. For systems such as in the scC0 2 phase with limited amount of water, the total carbonation 
potential is limited by the amount of water present for hydration of basalt. 



Introduction 

Underground sequestration of carbon dioxide is a poten- 
tially viable greenhouse gas mitigation option as it 
reduces the release rate of C0 2 to the atmosphere 
[1]. C0 2 can be trapped subsurface by four storage 
mechanisms: (1) structural and stratigraphic trapping; 
(2) residual C0 2 trapping; (3) solubility trapping; and 
(4) mineral trapping [2]. Mineral trapping has been con- 
sidered as the safest mechanism in long-term storage of 
C0 2 [3]. 
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Mineral storage of C0 2 in basaltic rocks is favored 
over siliciclastic reservoirs both by the higher abundance 
of divalent metal ions in basalt and the faster reactivity 
of basaltic glass or crystalline basalt [4]. Moreover, 
basalts such as the Columbia River flood basalts (CRB) 
are abundant and in many places close to C0 2 point 
source emissions [5]. During the last decade several 
flood basalts around the world have been mapped for 
the possibility of C0 2 storage, and possible candidates 
such as CRB in USA and the Deccan traps in India have 
been identified [4-6]. 

To be a candidate for C0 2 storage, the flood basalt 
must have a proper sealing and sufficient injectivity, the 
latter limited by the available connected pore space. In 
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flood basalts, the connected pore space is typically found 
at zones containing abundant vesicles or in breccias be- 
tween basalt flows. Because central zones of flood basalts 
commonly are dense and impermeable without vesicles, 
and flows are laterally continuous over large areas and 
commonly stacked vertically for hundreds of meters, 
flow units can act as seals [5]. The non-porous inner 
parts of flows may however be penetrated by networks 
of vertical fractures. These fractures can be open and 
conductive, or closed by mineralization and non- 
conductive. 

The main objectives of this study were to performe 
batch- and ID diffusion-reaction numerical simulations 
to determine the geochemical potential for secondary 
carbonate formation and to estimate the volume changes 
and the possibility of self-sealing following the basalt- 
C0 2 interactions. The CRB system was used as an ex- 
ample case and our results were compared to earlier 
reported laboratory experiments and numerical simula- 
tions of C0 2 -basalt interactions. As C0 2 stored under- 
ground will distribute spatially in the reservoir to give a 
range of reactive conditions, such as the potential of 
reactions by H 2 0 dissolved in supercritical C0 2 [7,8] or 
reactions in the H 2 0-rich phase from residually trapped 
C0 2 , we divided the simulations into three systems 
representing different parts of C0 2 storage: (1) basalt al- 
teration in the H 2 0-rich phase at constant C0 2 pres- 
sure; (2) basalt alteration in a H 2 0 saturated C0 2 phase, 
and (3) reactions at the boundary of the C0 2 plume 
where C0 2 diffuses into the aquifer from the boundary 
of the C02 plume (Figure 1). 

Methods 

All thermodynamic and kinetic calculations were per- 
formed using the geochemical code PHREEQC-2 [9]. This 
code is capable of simulating complex interactions be- 
tween dissolved gases, aqueous solutions, and mineral 




S3. On the boundary zone 

Figure 1 Sketch of possible reaction settings during C0 2 

storage in basalt. System 1 (SI) is close to the injector and 

contains a wet C0 2 (0.5 mole% H 2 0 at 100 bar and 40 C) with no 

residual water; system 2 (S2) is fully in the H 2 0 rich phase with C0 2 

diffusing in from the plume boundary; and system 3 (S3) is at the 

boundary of the C0 2 plume with both sufficient non-wetting C0 2 at 

a constant C0 2 partial pressure of 100 bar and with sufficient water 

wetting the mineral surfaces and available for reactions. 
\ J 



assemblages in batch and ID advection- diffusion-reaction 
mode. As the code can only model fully saturated systems, 
natural systems must be simplified to end-member situa- 
tions, such as given by constant pressure boundary condi- 
tions as may be the case close to underground C0 2 
plumes, or the assumption of packages (batches) of water 
reacting along a reaction path with a homogenous sedi- 
ment or rock body. Based on these limitations we divided 
the simulations into three systems representing different 
parts of C0 2 storage: (1) basalt alteration in the H 2 0-rich 
phase at constant C0 2 pressure; (2) basalt alteration in a 
H 2 0 saturated scC0 2 phase, and (3) reactions at the 
boundary of the C0 2 plume where C0 2 diffuses into the 
aquifer from the boundary of the scC0 2 plume (Figure 1). 
In the second case, we assumed that the C0 2 phase had 
swept through the systems and dried out residual water, 
giving only dissolved water in the scC0 2 phase. In this 
case an upper limit of carbonation potential was estimated 
as reactions were allowed to occur until (nearly) all water 
was consumed, passing the upper 2 mol/Kgw theoretical 
limit for the Truesdell- Jones activity model [9]. 

The standard state adopted in this study for the thermo- 
dynamic calculations was that of unit activity for pure 
minerals and H 2 0 at any temperature and pressure. For 
aqueous species other than H 2 0, the standard state was 
unit activity of the species in a hypothetical 1 molal solu- 
tion referenced to infinite dilution at any temperature and 
pressure. For gases, the standard state was for unit fugacity 
of a hypothetical ideal gas at 1 bar of pressure. All simula- 
tions used the llnl.dat database based on the thermo.com. 
V8.R6.230 dataset prepared at the Lawrence Livermore 
National Laboratory, with additions of thermodynamic 
data for those phases not present (see description below). 

C0 2 fugacity coefficients were estimated according to 
the modified Redlich-Kwong (SRK) equation of state 
[10] and the solubility was adjusted for by a poynting 
correction term {exp(v C o2(P sat - P)/RT) where v denotes 
molar volume, P pressure, R the universal gas constant 
and T absolute temperature) [11]. The density of C0 2 at 
40 C and 100 bar was approximated from Bachu and 
Stewart [12] to be 600 Kg/m 3 and the solubility of water 
in scC0 2 at the same conditions was approximated to 
0.5 mole% [13,14] 

The simulations were divided into batch simulations 
of the H 2 0 rich and C0 2 rich phases respectively, 
and ID diffusion of C0 2 in the H 2 0 rich phase to 
obtain information on the C0 2 -basalt interactions 
over a continuous range of C0 2 pressures. The latter 
was solved by PHREEQC using d t C — Di&^C + q , 
where C denotes molal (mol/Kgw) concentration, q 
denotes the sink term, subscripts t and x refer to 
derivatives in time and x-direction respectively, and 
an efficient diffusion coefficient D L of 0.45xl0" 9 m 2 /s 
was used for C0 2 [15] and all solutes. 
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Dissolution rates of minerals in the basaltic rock were 
calculated according to a kinetic equation taking into ac- 
count pH and the distance from equilibrium: 



(1) 



where S is the reactive surface area (m 2 ), k t are rate con- 
stants (moles/m 2 s), a H is the H + activity, n is the reac- 
tion order with respect to H + and OH", and £1 is the 
saturation state given by: 



n 



(AG r 



(2) 



Where AG r is the Gibbs free energy of the reaction, R 
is the gas constant, and T is absolute temperature. Reac- 
tion rate constants for crystalline basalt (pyroxenes and 
plagioclase) were obtained from Palandri and Kharaka 
[16], and pH dependencies were taken from the same 
source. The dissolution rate of basaltic glass was calcu- 
lated according to the expression suggested by Gislason 
and Oelkers [17]: 



^'(w) s (iF (i - n) 



(31 



where k + is the far-from-equilibrium dissolution rate co- 
efficient. The saturation state term 1-il was approxi- 
mated to 1 (i.e., rate independent to distance from 
equilibrium) supported by earlier numerical estimates of 
glass-C0 2 reactivity suggesting an approximately linear 
relation between time and reaction progress for basaltic 
glass [18]. This expression takes into account the effect 
of pH as well as the effect of the concentration of solutes 
such as fluoride as they complex with Al 3+ and reduce 
the Al 3+ activity [19]. The specific surface area for basalt 
(m 2 /g) was estimated by: 



0 



J sp 



(t>Ps V P 



(4) 



where the ratio A p /V p denotes the ratio between pore 
surface and pore volume (m 1 ), 0 is connected porosity, 



and p s (g/m 3 ) is the density of the basalt solid estimated 
from the fraction of the individual basalt components. A 
S sp value of 1.52 x 10" 5 m 2 /g basalt (= 0.137 m 2 /Kg 
water) was obtained for the CRB using an average basalt 
solid density of 2.93 x 10 6 g/m 3 with 10% connected pore 
space and a A p /V p ratio of 400 m" 1 [5]. The reactive sur- 
face area was calculated from the mass of the glass and 
minerals present according to: 



Si = MiVLiSspXr 



(5) 



where M and n are molar mass and moles of mineral i, 
and X r is the fraction of the total mineral surface that is 
reactive. As X r is highly uncertain and is suggested to 
vary by orders of magnitude [20,21], we used a value of 
0.1 for the base case and varied X r from 1 to 10" 3 . The 
use of mass or mass fractions of the individual basalt 
components to estimate the release rates of elements 
from the basalt is supported by a recent experimental 
study which suggests that release rates estimated from 
the sum of volume fractions of the constituent minerals 
are within one order of magnitude from measured values 
[22]. A list of kinetic parameters is given in Table 1. All 
secondary phases were allowed to form according to the 
local equilibrium assumption [23]. 

Changes in solid-phase volumes and porosities 0 caused 
by the mineral reactions were calculated according to: 



Vt 



total 



t=0 



(6) 



where (p t=0 is the initial porosity, n and v are moles and 
molar volume of mineral i respectively, and V to tai is the 
total volume of the system. 

The basalt was defined to consist of a mixture of glass 
and crystalline basalt with mineral and glass fractions 
chosen based on reported data from CRB [6,24,25]. To 
represent the crystalline basalt, plagioclase 
(Cao.5Nao.5AlL5Si2.5Og) and the pyroxenes augite 
(Ca 0 . 7 Feo.6 Mg 0 . 7 Si 2 O 6 ) and pigeonite (Cai.14Feo.64 
Mg 0 .22Si2O 6 ) were chosen. The hydrolysis equilibrium 
constants of these phases were estimated using the 
PHREEQC program assuming ideal solid solutions of 
the end-members enstatite, ferrosilite and wollastonite 
for the pyroxenes, and albite and anorthite for the 



Table 1 Kinetic parameters for dissolution of primary minerals based on empirical data given in Palandri and Kharaka 
[16] and for basaltic glass from [17] 





k H (mol/m 2 s) 


Ea H kJ/mol 


n H 


k N (mol/m 2 s) 


Ea N kJ/mol 


k oh (mol/m 2 s) 


Ea OH kl/mol n OH 


References 


Augite 


1 .58e-7 


78 


0.7 


1.07e-12 


78 






[16] 


Pigeonite 


1 .58e-7 


78 


0.7 


1.07e-12 


78 






[16] 


Feldspar 


1 .58e-9 


53.5 


0.541 


3.39e-12 


57.4 


4.78e-15 


59 -0.57 


[16] 


glass 


1e-10 


25.5 


1 










[17] 


Magnetite 


2.57e-9 


18.6 


0.279 


1.66e-11 


18.6 






[16] 
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plagioclase. Equilibrium constants for the solid solutions 
for temperatures up to 100 C were estimated with 
PHREEQC and from these data coefficients a to e for 
the PHREEQC built-in analytical expression (log K = a + 
bT + c/T + dlog w (T) + e/T 2 ) were estimated using non- 
linear regression in MATLAB. 

The glass composition (Cao.015Feo.095 Mgo.oe5Nao.025 
K0.01AI0.105 So.003Sio.5O1.35) was taken from [6] and 
modified by adding a small fraction of sulfur which is a 
common minor constituent of the CR basaltic glass [26] . 

The secondary mineral assemblage was chosen based 
on reports on basalt weathering [27-30], with additional 
carbonates that could potentially form at elevated C0 2 
pressures from the release of Fe, Mg and Ca. The anker- 
ite composition chosen for this work was CaFe 0 .6 Mg 0 . 4 
(C0 3 ) 2 which corresponds to a solid solution of 0.6 
ankerite (CaFe(C0 3 ) 2 ) and 0.4 dolomite (CaMg(C0 3 ) 2 ). 
Because ankerite (CaFe 0 .6 Mg 0 .4(CO 3 ) 2 ) was not listed in 
the thermodynamic database, we estimated values using 
the same approach as in [31]. The full list of secondary 
minerals is given in Table 2. 

To simulate the CRB-C0 2 interaction we used the 
average concentrations of solutes reported for the Grand 
Ronde Formation (Table 3). As supercritical C0 2 (scC0 2 
at T>31.1 C; P> 73.9 bar) is the preferred choice for 
C0 2 storage, based on higher density compared to gas- 
eous C0 2 , we simulated aqueous-phase basalt-C0 2 
interaction at a depth of 800 meters at a C0 2 pressure 
of 100 bar and temperatures of 40 to 100 C. The reactiv- 
ity of basalt and a H 2 0 saturated scC0 2 phase was 
simulated using an estimated 0.5 mol% H 2 0 and a C0 2 
density of 600 g/cc giving an initial mass of 0.003 Kg 
H 2 0 per 1 liter pore space. 

Results 

System 1: Basalt alteration in the H 2 0-rich phase at 
constant C0 2 pressure 

i) CRB mineral and glass dissolution and formation of 
secondary minerals 

Following the injection of C0 2 into the system, pH im- 
mediately decreased from 9.5 to below 4, and thereafter 
gradually increased to 5.8 at the end of 10000 years 
(Figure 2a). At the acidic pH secondary phases such as 
saponite (Ca 0 .i65Mg 3 Alo. 33 Si 3 . 6 70 10 (OH) 2 ), celadonite 
(KMgAlSi 4 O 10 (OH) 2 ) and zeolite (stilbite) were thermo- 
dynamically stable and formed (Figure 2b). Glass dis- 
solved orders of magnitude faster than the crystalline 
basaltic constituents and more than half dissolved after 
10000 years (Figure 2c). The dissolution rate of glass 
was not increased by the aqueous fluoride as the Al 3+ 
activity was fixed by the kaolinite and amorphous silica 
equilibria. The fluoride therefore only increased the total 
soluble aluminium. The steady release of Fe from the 
basalt saturated the water with respect to siderite and a 



total amount of 10 moles/kgw formed after 10000 years 
(Figure 2d). Other carbonates, such as ankerite, dolo- 
mite, magnesite, and dawsonite did not form as elements 
such as Mg and Ca was consumed by the non-carbonate 
secondary phases. 

The effect of temperature on the basalt hydration and 
carbonation was investigated by simulating the system at 
60, 80 and 100 C (Figure 3). As for the 40 C simulation 
we see that basaltic glass dissolves orders of magnitude 
faster than the crystalline basalt components and the 
glass is the major source for the secondary phases. The 
dissolution rates of the basalt components scale expo- 
nentially with temperature, and the glass is almost 



Table 2 Mineralogy included in the model 



Primary minerals 

1 Augite (En0.35Fs0.3Wo0.35) 
1 Pigeonite (En0.57Fs0.32Wo0.1 1) 
1 Plagioclase (An50) 
Glass Ca 0 .oi5Fe 0 .o95Mgo .065 

N3o.025Ko.01 AI0.1 05 So.OO3Sio.5Ol .35 

Magnetite 

Secondary minerals 

Si0 2 (am) 

Albite 

Goethite 

Calcite 

Hematite 

Kaolinite 

Smec high Fe-Mg 

Saponite-Mg 

Celadonite 

Stilbite 

Dawnsonite 

Siderite 

1 Ankerite (Ank 0 . 6 Do 0 .4) 

Dolomite 

Magnesite 



Initial Density 2,3 Log K° 

Weight % (g/cm 3 ) 



16 3.40 21.00 

3 3.38 21.40 

35 2.68 14.20 

45 2.92 -99.00 

1 5.15 10.47 

0 2.62 -2.71 

0 2.62 2.76 

0 3.80 0.53 

0 2.71 1.85 

0 5.30 0.11 

0 2.60 6.81 

0 2.70 1 7.42 

0 2.40 26.25 

0 3.00 7.46 

0 2.15 1.05 

0 2.42 4.35 

0 3.96 -0.19 

0 3.05 -19.51 

0 2.84 4.06 

0 3.00 2.29 



The mineralogy of the CRB has been described in [25,32] and the weight 
fraction of pyroxene, feldspar and glass was estimated as average values from 
the reported data. 

1 Solid solutions. En (enstatite), Fs (ferrosilite), Wo (wollastonite), An (anorthite), 
Ank (ankerite), Do (dolomite). For details on the calculations of the ankerite 
solid solution see [31]. 

Superscript 0 denotes Standard state (T = 298K, P = 1 atm). The equilibrium 
constant log K value is that for the forward dissolution reaction for one mole 
unit of the mineral. 

3 AII thermodynamic data {log K and coefficients for the PHREEQC analytical 
temperature expression) from the llnl.dat PHREEQC database, except for the 
solid-solutions first estimated in PHREEQC by ideal solid solutions and then 
added to the PHREEQC database as new solid solution phases. 
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Table 3 Composition of initial formation water 



Elements (totals) 


Mol/kgw 


Na 


1.0 x 10~ 3 


Ca 


6.0 x 10~ 4 


K 


1.0 x 10~ 4 


Mg 


2.0 x 10~ 5 


Fe 


1 .2 x 1 0" 6 


Alkalinity (HC0 3 ) 


2.0 x 10~ 3 


CI 


3.0 x 10~ 4 


S (SOf) 


1.0 x 10~ 4 


Si 


2.0 x 10~ 4 


Al 


1.0 x 10~ 6 


Log(0 2 ) 


-10.68 


PH 


7.5 



completely dissolved after 10000 years at 60 C, whereas 
the time for a complete dissolution takes 4000 and 
1500 years at 80 and 100 C respectively (Figure 3a, d, g). 
The secondary mineral assemblages were largely the 
same for all temperatures. Stilbite dominated together 
with amorphous silica (40 and 60 C) and quartz (80 and 
100 C) (Figure 3b, e, h). Saponite formed at 40 and 
60 C, but not at higher temperatures. Other secondary 
minerals such as albite, celadonite, and kaolinite formed 
at all conditions. At 60 C, magnesite and dolomite were 
still considered to be too slow to form (see [29]) and 



siderite was the only phase that formed. At 80 and 
100 C, magnesite and later ankerite formed together 
with siderite. Taking zero porosity as the maximum ex- 
tent of possible reactions we see that the total amount of 
C0 2 trapped as solid carbonates did not change much 
with temperature (Figure 4). The reaction rates 
increased however with temperature and the time 
needed to reach the maximum potential therefore 
decreased with temperature (Figure 4). 

ii) On the limitation of pore-space for the basalt 
carbonation 

Secondary phases such as stilbite and amorphous silica 
have lower density than the basalt components and al- 
teration therefore leads to a reduction of pore space. At 
the presence of C0 2 , secondary carbonates further re- 
duce the pore space. For the volume calculations we 
used expression (6) with the molar volumes listed in 
Table 2. At 40 C, the starting porosity of 10% is reduced 
to 0.85% after 10000 years. At the higher temperatures 
all porosity is lost after 2700, 1200, and 300 years re- 
spectively at 60, 80, and 100 C (Figure 5). Taking the ex- 
treme of 0% porosity as the limit for the reactions we 
obtain a maximum carbonation potential (mol C0 2 
stored/Kgw) at the different temperatures of 13.5, 29.3, 
and 28.5 moles for 60, 80, and 100 C (Figure 5). The 
simulated clogging of the pore space fits well with short- 
term laboratory percolation experiments on open-system 



a 



8 
7.5 

7 
6.5 

6 

; 5.5 
l 

5 
4.5 
4 

3.5 
3 




4000 6000 
Time, years 



8000 10000 



— Hematite 


— Clinoptilolite-Ca ^ 


Albite 


Stilbite 


Saponite 


Pyrite y^' 


Kaolinite 




Si02 




Celadonite 





g * 



14 
12 
10 
8 
6 
4 
2 
0 



2000 4000 6000 8000 10000 



2000 



8000 10000 



2000 



4000 6000 
Time, years 



8000 10000 



4000 6000 
Time, years 

Figure 2 Basalt alteration at 40 C and 100 bar C0 2 pressure over 10000 years, a) pH changes; b) mass fractions of basaltic glass and 
crystalline basalt components; c) secondary phases formed; and d) moles of secondary carbonates (siderite) formed per kgw. 
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glass; 
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Mass fractions of minerals following basalt alteration at 60, 80, and 100 C over 10000 years: a, d, g) primary basalt minerals and 
h) secondary phases except the carbonates; and c, f, j) carbonates. 



basalt-C0 2 alteration which shows loss of porosity and a 
rapid reduction of permeability during C0 2 -basalt inter- 
actions (e.g., [33]). 

Hi) Reduction of pore-space as a function of reactive surface 
area 

As the reactive surface area is a large uncertainty we 
simulated the changes of porosity over a range of values 
from a maximum being equal to the estimated physical 
surface area S 0 (equation (5) with X r = l) to a three 
orders of magnitude reduction (Figure 6). The physical 
conditions of the simulated system was the same as for 
the base-case at 40 C and a C0 2 pressure of 100 bars. 
At a reactive surface area that is equal to the estimated 



S 0 all porosity is lost after approximately 1000 years as 
stilbite and siderite fills the pore space. If the reactive 
surface area is reduced by one order of magnitude (i.e. 
the base case) nearly 1/10 of the original 10% porosity is 
preserved. Further reductions by one and two orders of 
magnitude lead to smaller changes and at three orders 
of magnitude reduction relative to S 0 almost no change 
is observed (Figure 6). 

System 2: The potential for carbonate growth in a H 2 0- 
saurated scC0 2 phase 

The reaction between H 2 0 dissolved in scC0 2 and bas- 
alt was simulated at 100 bar pressure and 40 C. The ini- 
tial amount of water was 0.003 Kg and no H 2 0 was 
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1 100 10000 

Time, years 

Figure 4 C0 2 trapped in solids for 40 to 100 C simulations at 
100 bar C0 2 pressure. A cut-off value is used when all pore space 
is filled up with the secondary phases (see Figure 5). The simulations 
suggest that the total amount of secondary carbonates that form is 
dictated by the available pore space and the thermodynamic 
stability of secondary phases rather than temperature, whereas 
carbonate generation rates depend on the exponential increase of 
basalt dissolution rates with temperature. 



allowed to enter the system. This is an ideal end- 
member case and serves to illustrate the carbonation po- 
tential in a volume with limited hydration potential. 

As secondary phases such as stilbite formed, water 
was rapidly consumed and most gone after 45 years 




0 2000 4000 6000 8000 10000 
Time, years 



Figure 5 Porosity changes caused by the basalt alteration at 40 

to 100 C. Secondary hydrated species and carbonates incorporate 
the H 2 0 and C0 2 masses into the solids and clogs the pore space. 
As reaction rates increase exponentially with temperature, the pore 
space is filled up faster at the higher temperatures. 
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Figure 6 Porosity changes caused by the basalt alteration at 

40 C and specific surface areas ranging from S Q (estimated 

total basalt surface area) down to a three orders of magnitude 

reduction. The base case specific surface used was So/10, 
k ) 

(Figure 7a). At this point stilbite was unstable and sup- 
plied water until all water was consumed after approxi- 
mately 100 years (Figure 7a). Following the basalt 
hydration, siderite and ankerite formed from the 
released Ca, Mg, and Fe, with a final total amount of 0.2 
moles C0 2 consumed per liter pore space after 100 years 
(Figure 7b). If H 2 0 had been allowed to dissolve into the 
scC0 2 phase from residual aqueous phases trapped in 
the smaller pores, the carbonation potential would have 
been larger. This process is however, to the knowledge 
of the authors, not possible to simulate using the 
PHREEQC code, and was therefore outside the scope of 
this study. 

System 3: ID diffusion of C0 2 into the CRB aquifer 

To see how the basalt reacted under different C0 2 pres- 
sures, we defined a ID diffusion-reaction simulation. 
This provided us with basalt-C0 2 interactions over a 
continuous range of C0 2 pressures from the background 
1 bar up to the maximum 100 bars. The system corre- 
sponds to a stagnant zone presented as a column with 
one end close to the boundary of the injected C0 2 
plume and the other end further away from the plume 
(Figure 1). The distance reached for the C0 2 into the 
column is given by the balance between diffusive trans- 
port and removal of carbon by secondary carbonate for- 
mation. We therefore varied reaction rates from no 
reaction giving the maximum lengtht of diffusive trans- 
port, and up to the base-case rate given by a reaction 
surface area 1 order of magnitude lower than the esti- 
mated physical surface area. Figure 8 shows pH, dis- 
solved C0 2 (mol/Kgw) and amount of secondary 
carbonate formed in the ID column. As C0 2 diffuses 
into the column pH drops to approximately 4 at full 
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saturation. The depth of diffusion into the ID column is 
approximately 40 meters at 1000 years if no carbonate 
forming reactions occur (Figure 8a). The penetration 
depth decreased rapidly if reactions were allowed as sid- 
erite formed and pulled carbon out of aqueous system 
(Figure 8b, c). At the highest reaction rate (base-case), 
C0 2 diffused less than 10 meters into the column as ap- 
proximately 13 moles/Kgw of siderite formed at the end 
of the 10000 years simulation (Figure 8d). Siderite 
formed at greater depth if the reactive surface area was 
reduced by another order of magnitude, but less formed 
in total (Figure 8d). 

Discussion 

Uncertainty on the reactive surface area 

The reactive surface area is considered as a major source 
of uncertainty (e.g., [20,34]) and this leads to corre- 
sponding high uncertainties in timing and extent of 
reactions as dissolution rates have a first order depend- 
ence on reactive surface areas. Weathering rates in na- 
ture are commonly observed to be 1-3 orders of 
magnitude lower than in laboratory experiments (e.g., 
[20,21,34]), and this may in part be explained by differ- 
ences in reactive and physical (total) surface area be- 
tween experimental and natural systems. We assumed in 
this study a base-case reactive surface area 1 order of 
magnitude lower than the estimated physical surface 
area for the basalt. A further two orders of magnitude 
reduction in the reactive surface area, which is within 
the range of values expected for natural systems, 
resulted in little basalt alteration and only minor reduc- 
tion of porosity (see Figure 5). A better understanding of 
the surface area of porous basalt and the effect of time 
(aging) on features such as dislocation densities and re- 
active surface areas are therefore required to understand 
the potential for C0 2 mineral storage in basaltic rocks. 



Uncertainty on the choice of secondary phases used in 
the model 

Growth rate experiments of carbonates such as magne- 
site and dolomite have shown that the activation energy 
is high and that growth is negligible at low temperatures 
(e.g., [35-37]). Dissolution rate studies of siderite sug- 
gests that the reaction rate is intermediate between cal- 
cite and magnesite [38,39], and growth rate data suggest 
that siderite may form down to room temperature [40]. 
Data on ankerite dissolution and growth is to the know- 
ledge of the authors not known. The crystallographic 
and physical characteristics of ankerite do resemble 
those of dolomite and siderite, and the chemistry is 
related to dolomite with the Mg 2+ substituted by various 
amounts of Fe 2+ and Mn 2+ . If the growth rate is close to 
the magnesian carbonates such as dolomite and magne- 
site [41,42], the amount that may form during low- 
temperature alteration is likely low. In this case, more 
iron would be available for siderite growth. If on the 
other hand the growth rate is closer to siderite, we 
would expect ankerite or other FeMg solid solution car- 
bonates to grow during low-temperature alteration. 

One uncertainty related to the local- equilibrium as- 
sumption is on the growth retention time for the sec- 
ondary carbonates. The local-equilibrium assumption 
predicts growth of the secondary phases as soon as an 
infinitesimally small supersaturation is reached [23]. The 
time it takes to nucleate sufficient mass to initiate a sig- 
nificant growth may however be hundreds to thousands 
of years for some secondary phases [31]. There are no 
nucleation rate data for siderite and ankerite and the re- 
tention time is hence unknown. 

Finally, the total potential for secondary carbonate 
growth may be affected by the amount of magnesium 
and iron that enters ferromagnesian calcites. As a signifi- 
cant fraction of the metal cations may substitute for cal- 
cium (e.g., [43]), a iron-magnesium rich calcite may 
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Figure 8 1 D reaction-diffusion of C0 2 into permeable basalt. 

Partial pressure of the inlet boundary was fixed at 100 bar and with 
a column temperature of 40 C. As consumption of C0 2 by siderite 
growth affects the depth of C0 2 diffusion, we ran a sensitivity study 
on reactive surface area going from no reaction (a) up to the base 
case (c). Finally the amount of C0 2 trapped as solid carbonate 
(siderite) was compared for the base case and reduced reactive 
surface area (d). We see that the reaction rates strongly constraint 
the depth of the column affected by the C0 2 diffusion. 



potentially form rather than ankerite and thereby reduce 
the amount of siderite formed. 



Comparisons to experiments, numerical simulations and 
natural analogues of basalt-C0 2 interactions 

Our simulations suggest that the potential for carbonate 
growth is limited to siderite or FeMg carbonates at low 
temperatures as secondary phases such as zeolites out- 
competed the carbonates for calcium. We here compare 
our simulated results with reported data on C0 2 -basalt 
interactions from laboratory experiments, natural analo- 
gues, and other reported numerical simulations. 

The reactivity of CRB and other continental flood 
basalts are available from the long-term (months to 
years) laboratory experiments done by Schaef and co- 
workers [6,24]. In these experiments basalt samples from 
USA, India, South Africa, and Canada were reacted with 
C0 2 at about 100 bars and 60 to 100 C. Reacted samples 
from these experiments showed generation of Ca-rich 
carbonates interpreted as calcites with minor siderite 
and magnesite. In experiments on CRB using mixtures 
of H 2 S and C0 2 at 60 C and 100 bar and run for 
181 days, pyrite (FeS 2 ) formed together with Mg-Fe poor 
calcite and a Ca-poor Fe-carbonate [6]. Our simulations 
at the same temperatures show rapid formation of sider- 
ite (60 C) or siderite and magnesite at higher tempera- 
tures (Figure 3). Our simulations do not predict any 
calcite growth as the calcium activity is lowered by zeo- 
lite formation. Calcite would however form in our mod- 
els if the zeolites were not allowed to form at local 
equilibrium, and possibly if a magnesian ferroan (solid 
solution) calcite was used in the model instead of the 
pure end-member calcite. Therefore, the apparent differ- 
ence between our model and the experiment may be 
caused by our use of the local equilibrium assumption, 
whereas the zeolites in the laboratory experiments did 
not form at low temperatures due to slow kinetics. Re- 
cent experiments on basalt dissolution support the pre- 
ferential release of Mg and Fe over Ca at acidic 
conditions [22], suggesting that the MgFe-carbonates 
will dominate as secondary carbonates during C0 2 stor- 
age in basalt. 

Our numerical simulations share some similarities to 
other works such as by Marini [18] and Gysi [44], but 
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our model and hence the outcome is different in several 
aspects. The most comprehensive work done earlier is 
the numerical simulations done by Marini [18] on the 
reactivity of crystalline and glassy CFB following C0 2 
storage. The initial mineralogy was similar to our study 
whereas the temperature of 60 C was slightly higher 
than our base case 40 C. In [18] the C0 2 -basalt interac- 
tions were stretched to last for more than 280000 years 
compared to our 10000 years perspective. The main dif- 
ferences between our model and [18] are on the choice 
of secondary mineral assemblage, and on the focus of 
limiting factors such as the availability of water for hy- 
dration in the present work. The lack of zeolites and hy- 
drous phases other than kaolinite and goethite in [18] 
made Ca available for secondary carbonates and the total 
potential for carbonate formation was higher than in our 
work. Marini allowed dolomite and magnesite to form at 
60 C, whereas our simulations only produced siderite at 
the similar conditions. Moreover, the formation of daw- 
sonite in [18] is still uncertain and possibly limited at 
high silica activities and with an assemblage of stable 
NaAl-silicates defined to form [45]. Based on two differ- 
ent approaches, the reactive surface area for basalt 
was estimated to quite similar values. We estimated a 
specific surface area of approximately 1.5xl0~ 5 m 2 /g basa i t 
(= 0.14 m 2 /Kg water at 10% porosity) based on the Ap/V p 
values estimated by [46] and reported in [5], and reduced 
this value by one order of magnitude to get the reactive 
surface area. Marini used a geometric model giving a 
reactive surface area of 0.41 m 2 /Kg water. The higher 
reactive surface area and higher temperature of [18] 
resulted in faster reactions and more rapid clogging of the 
pore space (within a few years). Studies of natural basalt 
systems at similar or higher temperatures may give some 
insight into how fast pore space is clogged by basalt 
hydration or carbonation, and this should be used to 
improve the estimates of reactive surface areas of basalt 
for future studies. 

Another numerical study on low-temperature (25 C, 
30 bar C0 2 ) basaltic glass alteration was presented by 
Gysi et al. [44]. Again a main difference is on the choice 
of secondary minerals. Gysi et al. [44] allowed dolomite, 
magnesite, and Fe-Mg carbonate to form together with 
calcite and siderite, whereas we did not allow other Mg- 
Fe carbonates to form than ankerite. As previously sta- 
ted, the low-temperature formation of dolomite and 
magnesite is not likely because of the high apparent acti- 
vation energy and small kinetic coefficients for the 
growth of Mg-carbonates [35-37]. Other carbonates such 
as siderite and potentially FeMg-calcites are more likely 
to form at these low temperatures. The high reactive 
surface area used in [44] is based on a geometric model 
for glass fragments, and is hence not directly comparable 
with the surface area estimated for a vesicle pore space 



of a solid basalt. Although no inverse modeling was 
done to estimate the reactive surface area of the basalt 
in [44], fragmented basaltic rocks such as hyaloclastite 
breccias are expected to have significantly higher react- 
ive surface areas than porous solid basalts, and they are 
therefore correspondingly more reactive. 

One example of a natural analogue that shed light on 
C0 2 basalt interactions is the C0 2 charged basalt hosted 
groundwaters at Hekla, Iceland. Solution aqueous spe- 
cies sampled from natural cold springs and rivers here 
showed a drop in total inorganic carbon (TIC) that was 
interpreted to result from considerable formation of sec- 
ondary carbonate phases such as calcite [47]. Reaction 
path modeling of the system suggests however that the 
carbonate formation is associated with high pH in ac- 
cordance with the low TIC in the sampled waters. This 
system is therefore different from basalt C0 2 storage 
projects where higher C0 2 pressures may be maintained 
over time and the pH is lower. In addition to calcite, 
dolomite was also suggested as a potential storage host 
for the low temperature reactions in Hekla [47]. This 
may however be questionable as long-term laboratory 
experiments at room temperature have failed to form 
dolomite even at significant super saturations [48], 
explained by the high activation energy for dolomite 
growth [32,41]. Another natural analogue that more 
closely corresponds to industrial C0 2 storage is the 
basalt-hosted petroleum reservoir on Nuussuaq, West 
Greenland. In this system the bulk carbonate formation 
appears to have occurred as secondary weathering pro- 
ducts. Other alteration products such as zeolites and 
oxides were replaced by dolomite, magnesite, siderite, 
and calcite at temperatures of 70-120 C [49]. Therefore, 
taking into account the basalt weathering products and 
not only primary basalt minerals appears to be vital in 
estimating the total potential for secondary carbonate 
formation and the long-term potential for C0 2 storage 
in basalt systems. 

Summary and conclusions 

Simulations of closed-system (Pco2 = 100 bar, 40 C) and 
ID reaction-diffusion (P CQ2 = 0-100 bar, 40 C) alter- 
ation of basalt suggest that the potential of secondary 
carbonate formation is limited to siderite at low tem- 
peratures as divalent metal cations are preferentially 
consumed by zeolites and oxides. Higher temperatures 
60 - 100 C appear to be in favor of secondary carbonate 
formation, allowing the precipitation of carbonates such 
as magnesite, siderite and possibly dolomite and other 
FeMg carbonates (ankerite). Given an unlimited source 
of C0 2 (fixed C0 2 pressure), the total amount of C0 2 
stored as solid carbonates is orders of magnitude higher 
than the 1-2 mol/Kg water solubility in the formation 
water (Figure 4). The total amount trapped might 
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however be reduced if C0 2 , H 2 0 or pore space are limit- 
ing factors. The formation of secondary hydrous and 
carbonate phases increases the volume of solids and the 
porosity is correspondingly reduced (Figure 5). This to- 
gether with the immobilization of C0 2 by solid carbon- 
ate formation is in favor of safe long-term storage of 
C0 2 in basaltic aquifers. 

Competing interests 

We have no competing interests with any organization to publish the 
manuscript 'On the potential for C02 mineral storage in continental flood 
basalts - PHREEQC batch- and 1D diffusion-reaction simulations'. 

Authors' contributions 

VTH Pham has made substantial contributions to conception and designs 
the manuscript. Her contributions are to acquisition of data, analysis and 
interpretation of data. PA did revising and has given final approval of the 
version to be published. HH involved in drafting the manuscript, writing 
methodology part and revising of the final version. All authors read and 
approved the final manuscript. 

Acknowledgements 

We highly appreciated constructive comments and suggestions from the 
reviewers. This work has been funded by SSC-Ramore (Subsurface storage of 
carbon dioxide - risk assessment, monitoring and remediation) project and 
(partially) by SUCCESS centre for C0 2 storage under grant 193825/S60 from 
Research Council of Norway (RCN). SUCCESS is a consortium with partners 
from industry and science, hosted by Christian Michelsen Research as. 

Received: 1 June 201 1 Accepted: 14 June 2012 
Published: 14 June 2012 

References 

1. Holloway S: Underground sequestration of carbon dioxide-a viable 
greenhouse gas mitigation option. Energy 2004, 30:2318-2333. 

2. Bachu S, Bonijoly D, Bradshaw J, Burruss R, Holloway S, Christensen NP, 
Mathiassen OM: C0 2 storage capacity estimation: Methodology and gaps. 

Int J Greenhouse Gas Control 2007, 1 :430-443. 

3. IPCC: IPCC Special Report on Carbon Dioxide Capture and Storage. In In 

Prepared by Working Group III of the Intergovernmental Panel on Climate 
Change. Edited by Metz B, Davidson 0, Coninck H, Loos M, Meyer L. 
Cambridge, United Kingdom and New York, NY, USA: Cambridge University 
Press; 2005:442 pp. 

4. Oelkers EH, Gislason SR, Matter J: Mineral Carbonation of C02. Elements 
2008, 4:333-337. 

5. McGrail BP, Schaef HT, Ho AM, Chien Y-J, Dooley JJ, Davidson CL: Potential 
for carbon dioxide sequestration in flood basalts. J Geophysical Res-Solid 
Earth 2006, 1 1 1:1-13. 

6. Schaef HT, McGrail BP, Owen AT: Carbonate mineralization of volcanic 
province basalts. Int J Greenhouse Gas Control 2010, 4:249-261. 

7. Schaef HT, Windish CF, McGrail BP, Martin PF, Rosso KM: Brucite [Mg(OH) 2 ] 
carbonation in wet supercritical C0 2 : An in situ high pressure x-ray 
diffration study. Geochim Cosmochim Acta 201 1, 75:7458-7471. 

8. White MD, McGrail BP, Schaef HT, Hu JZ, Hoyt DW, Felmy AR, Rosso KM, 
Wurstner SK: Multiphase sequestration geochemistry: Model for mineral 
carbonation. Energy Procedia 201 1, 4:5009-5016. 

9. Parkhurst DL, Appelo CAJ: User's guide to PHREEQC (version 2) - a computer 
program for speciation, reaction-path, ID-transport, and inverse geochemical 
calculations:. U.S. Geological Survey, Water-Resources Investigation Report; 
1999:312. 

10. Soave G: Equilibrium constants from a modified Redlich-Kwrong 
equation of state. Chem Eng Sci 1972, 27:1 197-1203. 

1 1 . Hellevang H, Kvamme B: ACCRETE - Geochemistry solver for C0 2 -water- 
rock interactions. Proceedings GHGT 8 conference 2006, :8p. 

12. Bachu S, Stewart S: Geological sequestration of anthropogenic carbon 
dioxide in the western Canada sedimentary basin: Suitability analysis. 
J Can Pet Technol 2002, 41:32-40. 



13. Coan CR, King AD: Solubility of Water in Compressed Carbon Dioxide, 
Nitrous Oxide, and Ethane - Evidence for Hydration of Carbon Dioxide 
and Nitrous Oxide in Gas Phase. J Am Chem Soc 1971, 93:1857-1862. 

14. Sabirzyanov AN, llin AP, Akhunov AR, Gumerov FM: Solubility of water in 
Supercritical carbon dioxide. High Temp 2002, 40:203-206. 

15. Bahar M, Lui K: Measurement of the diffusion coefficient of C0 2 in formation 
water under reservoir conditions: Implication for C0 2 storage, Society of 
Petroleum Engineers, SPE Asia Pacific Oil and Gas Conference and Exhibition, 
20-22 October 2008. Perth, Australia:; 2008. 

16. Palandri JL, Kharaka YK: A Complilation of Rate Parameters of Water-Mineral 
Interaction Kinetics for Applicatin to Geochemical Modeling. U.S. Geological 
survey, Open report.: ; 2004:1068. 

17. Gislason SR, Oelkers EH: Mechanism, rates, and consequences of basaltic 
glass dissolution: II. An experimental study of the dissolution rates of 
basaltic glass as a function of pH and temperature. Geochim Cosmochim 
Acta 2003, 67:3817-3832. 

18. Marini L: Chapter 7: Reaction path Modelling of geological C0 2 sequatration, 
In Geological Sequestration of Carbon Dioxide Thermodynamics, Kinetics, and 
Reaction Path Modeling:. Dev Geochem; 2006:319-394. 

1 9. Wolff-Boenisch D, Gislason SR, Oelkers EH: The effect of fluoride on the 
dissolution rates of natural glasses at pH 4 and 25 degrees C. Geochim 
Cosmochim Acta 2004, 68:4571-4582. 

20. White AF, Peterson ML: Role of reactive-surface-area characterization in 
geochemical kinetic models. In In. Chemical Modeling of Aqueous Systems II. 
Edited by Melchior D, et al. Washington, DC: ACS Symposium Series, ACS; 1990. 

21 . White AF, Brantley SL: The effect of time on the weathering of silicate 
minerals: why do weathering rates differ in the laboratory and field? 
Chem Geol 2003, 202:479-506. 

22. Gudbrandsson S, Wolff-Boenisch D, Gislason SR, Oelkers EH: An experimental 
study of crystalline basalt dissolution from 2 < pH < 1 1 and temperatures 
from 5 to 75 C. Geochim Cosmochim Acta 201 1, 75:5496-5509. 

23. Helgeson HC: Evaluation of Irreversible Reactions in Geochemical 
Processes Involving Minerals and Aqueous Solutions .I. Thermodynamic 
Relations. Geochim Cosmochim Acta 1 968, 32:853-877. 

24. Schaef HT, McGrail BP, Owen AT: Basalt-C0 2 -H 2 0 Interactions and Variability 
in Carbonate Mineralization Rates. Energy Procedia 2009, 1:4899-4906. 

25. Schaef HT, McGrail BP: Dissolution of Columbia River Basalt under mildly 
acidic conditions as a function of temperature: Experimental results 
relevant to the geological sequestration of carbon dioxide. Appl Geochem 
2009, 24:980-987. 

26. Blake S, Self S, Sharma K, Sephton S: Sulfur release from the Columbia 
River Basalts and other flood lava eruptions constrained by a model of 
sulfide saturation. Earth Planet Sci Lett 2010, 299:328-338. 

27. Neuhoff PS, Fridriksson T, Arnorsson S, Bird DK: Porosity evolution and 
mineral paragenesis during low-grade metamorphism of basaltic lavas 
at Teigarhorn, eastern Iceland. Am J Sci 1999, 299:467-501. 

28. Neuhoff PS, Rogers KL, Stannius LS, Bird DK, Pedersen AK: Regional very 
low-grade metamorphism of basaltic lavas, Disko-Nuussuaq region, West 
Greenland. Lithos 2006, 92:33-54. 

29. Stefansson A, Gislason SR, Arnorsson S: Dissolution of primary minerals in 
natural waters - II. Mineral saturation state. Chem Geol 2001, 172:251-276. 

30. Reidel SP, Johnson VG, Spane FA: Natural gas storage in basalt aquifers of 
the Colombia Basin, Pacific Northwest USA: A guide to site characterization. H: 
Report Pacific Northwest National Laboratory; PNNL-13962; 2002. 

31 . Pham VTH, Lu P, Aagaard P, Zhu C, Hellevang H: On the potential of C0 2 - 
water-rock interactions for C0 2 storage using a modified kinetic model. 
International Journal of Greenhouse Gas Control 201 1, 5:1002-1015. 

32. Reidel SP: A lava flow without a source: the cohassett flow and its 
compositional components, Sentinel Bluffs Member, Columbia River 
Basalt Group. J Geol 2005, 1 13:1-21. 

33. Peuble S, Godard M, Gouze P, Luquot L: C0 2 sequestration in olivine rich 
basaltic aquifers: a reactive percolation experimental study. Geochim 
Cosmochim Acta 2009, 73:A1635. 

34. Velbel MA: Constancy of silicate-mineral weathering-rate ratios between 
natural and experimental weathering: implications for hydrologic control 
of differences in absolute rates. Chem Geol 1993, 105:89-99. 

35. Saldi GD, Jordan G, Schott J, Oelkers EH: Magnesite growth rates as a 
function of temperature and saturation state. Geochim Cosmochim Acta 
2009, 73:5646-5657. 

36. Arvidson RS, Mackenzie FT: The dolomite problem; control of precipitation 
kinetics by temperature and saturation state. Am J Sci 1999, 299:257-288. 



Van Pham et al. Geochemical Transactions 2012, 13:5 
http://www.geochemicaltransactions.eom/content/13/1/5 



Page 12 of 12 



37. Arvidson RS, Mackenzie FT: Tentative kinetic model for dolomite 
precipitation rate and its application to dolomite distribution. Aquat 
Geochem 1996, 2:273-298. 

38. Golubev SV, Bebnezeth P, Schott J, Dandurand JL, Castillo A: Siderite 
dissolution kinetics in acidic aqueous solutions from 25 to 100 C and 0 
to 50 atm pC0 2 . Chem Geol 2009, 265:13-19. 

39. Pokrovsky OS, Golubev SV, Schott J, Castillo A: Calcite, dolomite and 
magnesite dissolution kinetics in aqueous solutions at acid to 
circumneutral pH, 25 to 150 C and 1 to 55 atm pC0 2 : New constraints 
on C0 2 sequestration in sedimentary basins. Chem Geol 2009, 265:20-32. 

40. Greenberg J, Tomson M: Precipitation and dissolution kinetics and 
equilibria of aqueous ferrous carbonate vs temperature. Appi Geochem 
1992, 7:185-190. 

41. Arvidson RS, Mackenzie FT: Tentative Kinetic Model for Dolomite 
Precipitation Rate and Its Application to Dolomite Distribution. Aquat 
Geochem 1997, 2:273-298. 

42. Saldi GD, Jordan G, Schott J, Oelkers EH: Magnesite growth rates as 
function of temperature and saturation state: An HAFM study. Geochim 
Cosmochim Acta 2009, 73:A1 149. 

43. Busenberg E, Plummer LN: Thermodynamics of magnesian calcite solid- 
solutions at 25 C and 1 atm total pressure. Geochim Cosmochim Acta 
1989,53:1189-1208. 

44. Gysi AP, Stefansson A: Numerical modelling of C0 2 -water-basalt 
interaction. Mineral Mag 2008, 72:55-59. 

45. Hellevang H, Declercq J, Aagaard P: Why is dawsonite absent in C0 2 
charged reservoirs? Oil & Gas Science and Technology - Re. IFP Energies 
nouvelles 201 1,66:1 19-135. 

46. Saar MO, Manga M: The relationship between permeability, porosity and 
microstructure in vesicular basalts.: Master Thesis, Univ. of Oregon, Eugene; 
1998:91. 

47. Flaathen TK, Gislason SR, Oelkers EH, Sveinbjornsdottir AE: Chemical 
evolution of the Mt. Hekla, Iceland, groundwaters: a natural analogue for 
C0 2 sequestration in basaltic rocks. Appi Geochem 2009, 24:463-474. 

48. Land LS: Failure to precipitate dolomite at 25 degrees C from dilute 
solution despite 1000-fold oversaturation after 32 years. Aquat Geochem 
1998, 4:361-368. 

49. Rogers KL, Neuhoff PS, Pedersen AK, Bird DK: C0 2 metasomatism in a 
basalt-hosted petroleum reservoir, Nuussuaq, West Greenland. Lithos 2006, 

92:55-82. 



doi:1 0.1 186/1467-4866-13-5 

Cite this article as: Van Pham et al.: On the potential for C0 2 mineral 
storage in continental flood basalts - PHREEQC batch- and 1D 
diffusion-reaction simulations. Geochemical Transactions 2012 13:5. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at (^\ RiftMM i rpntral 

www.biomedcentral.com/submit \^ ™omea centra I 



